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ABSTRACT 

We confirm the planetary nature of Kepler-539 b (aka Kepler object of interest K00372.01), a giant transiting exoplanet 
orbiting a solar-analogue G2 V star. The mass of Kepler-539 b was accurately derived thanks to a series of precise radial 
velocity measurements obtained with the CAFE spectrograph mounted on the CAHA 2.2-m telescope. A simultaneous 
fit of the radial-velocity data and Kepler photometry revealed that Kepler-539 b is a dense Jupiter-like planet with 
a mass of Mp = 0.97 ± 0.29 Mjup and a radius of Rp = 0.747 ± 0.018 J?jup, making a complete circular revolution 
around its parent star in 125.6 days. The semi-major axis of the orbit is roughly 0.5 an, implying that the planet 
is at « 0.45 au from the habitable zone. By analysing the mid-transit times of the 12 transit events of Kepler-539 b 
recorded by the Kepler spacecraft, we found a clear modulated transit time variation (TTV), which is attributable to 
the presence of a planet c in a wider orbit. The few timings available do not allow us to precisely estimate the properties 
of Kepler-539 c and our analysis suggests that it has a mass between 1.2 and 3.6Mjup, revolving on a very eccentric 
orbit (0.4 < e < 0.6) with a period larger than 1000 days. The high eccentricity of planet c is the probable cause of 
the TTV modulation of planet b. The analysis of the CAFE spectra revealed a relatively high photospheric lithium 
content, A(Li) = 2.48 ± 0.12 dex, which, together with both a gyrochronological and isochronal analysis, suggests that 
the parent star is relatively young. 
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1. Introduction 


The transiting extrasolar planet (TEP) population turns 
out to be the best atout in the hands of exoplanetary scien- 


(e-g-, 

Bakos et al.||2004 

Alonso et al.||2004 

McCullough et 

al. 1120051 IPollacco et ai. 

2006 

IPepper et al.||2007t |Alsubai 

et al. 

2013 Bakos et al. 

2013 

I and then to those from the 

space 

(CoRoT: Barge et ai. 

2008, Kepler-. 

Borucki et al. 


20111, more than 1290 transiting planetary systems have 


now been found. The possibility to easily derive most of 
their physical and orbital parameters, and to investigate 
even the composition of their atmosphere, make the TEPs 
the most suitable targets to detect and study in detail. They 
currently represent the best statistical sample to constrain 
the theoretical models of planetary formation and evolu¬ 
tion. Moreover, transit time variation (TTV) studies al¬ 
lowed the detection of additional non-transiting bodies in 
many TEP systems, highlighting the great efficacy of high- 


precision and high-cadence photometry in the search for 
extrasolar planets. 

The majority of TEPs have been discovered by the 
Kepler space telescope, which has revealed how varied they 
are, in terms of mass and size, and how diverse their ar¬ 
chitectures can be, confirming science-fiction pictures in 
several cases or going beyond human imagination in oth¬ 
ers. The numerous TEP discoveries achieved by Kepler are 
based on time-series photometric monitoring of the bright¬ 
ness of over 145 000 main sequence stars, in which the 
periodic-dimming signal caused by a TEP can be found. 
However, such a signal can be mimic by other astrophysical 
objects in particular configurations. Radial velocity (RV) 
follow-up observations of the ‘potential’ parent stars are 
therefore a fundamental step for distinguishing real TEPs 
from false positive cases. 

Here we focus our attention on the Kepler system 
Kepler-539. Thanks to precise RV measurements, presented 
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in Sect.[^ we confirm the planetary nature of Kepler-539 b 
(aka KOI-372b, K00372.01, KIC 6471021), a dense gas- 
giant planet moving on a wide orbit around a young and 
active G2V star {V = 12.6 mag). The analysis of the phys¬ 
ical parameters of Kepler-539 b and its parent star are de¬ 
scribed in Sect.[^ Moreover, the clear variation observed in 
the mid-times transits of Kepler-539 b strongly support the 
existence of an additional massive planet, Kepler-539 c, in 
the system, moving on a larger and very eccentric orbit, as 
discussed in Sect. HI We summarise our results in Sect.[5l 


2. Observations and data analysis 

2.1. Kepler photometry 

The Kepler spacecraft monitored Kepler-539 from quar¬ 
ters 0 to 17 (i.e. four years; from May 2009 to May 2013). 
It was labelled as a Kepler object of interest (KOI) due 
to a ~ 0.2% dimming in its light curve with a period 
of ^ 125days ( Borucki et al.| 2011|. This periodic dim¬ 
ming is actually caused by the transit of a Jupiter-like 
planet candidate, Kepler-539 b, moving on a quite wide or¬ 
bit around the star. 12 transits of Kepler-539 b are present 
in the Kepler long cadence (LC) light curves. We have 
labelled them from cycle -5 to cycle 6 (see Fig.[^. Two 
of the transits are incomplete (cycles -2 and -1); two are 
most likely contaminated by starspot crossing anomalies 
(cycles -5 and 2); three were also covered in short cadence 
(SC) (cycles 4, 5 and 6; see Fig.[^. The complete Kepler 
light curve is shown in Fig.[^ which highlights a signifi- 
cant stellar variability (0.0470 ± 0.0002 mag peak-to-peak). 
McQuillan et al. (20131 found a periodic photometric mod- 
ulation in the light curve and, by assuming that it is in¬ 
duced by a star-spot activity, estimated a stellar rotation 
period of 11.769 ± 0.016 days. This value is in good a gree- 
ment with that found by Walkowicz & Basri (2013a|b), i.e. 
11.90 ±3.45 days. 


2.2. Radial velocity follow-up observations 

We monitored Kepler-539 between July 2012 and July 
2015 wi th the Calar Alto Fib re-fed Echelle spectrograph 
(CAFE; jAceituno et al. 2013) mounted on the 2.2 m tele¬ 
scope at the Calar Alto Observatory (Almen'a, Spain) as 
part of our follow-up programme of Kepler candidates. 
This programme ha s already confirmed the planetary na- 
ture of Kepler-91 b (Lillo-Box et al.| 2014a|b|), Kepler-432 b 


(Ciceri et al. 2015), Kepler-447b (Lillo-Box et al. 2015bI 


and has identihed and characterised some falsepositives 


in the sa mple of Kepler planet candidates (Lillo-Box et 
al. 2015a). CAFE has a spectral coverage from 4000 A to 
9500 A, divided into 84 orders with a mean resolving power 
of i? = 63 000. We acquired 28 spectra, which were re¬ 
duced usi ng the dedicated pipe line provided by the ob¬ 
servatory (Aceituno et al. 2013). We used thorium-argon 
(ThAr) exposures obtained after each science spectrum 
to wavelength-calibrate the corresponding data. The final 
spectra have signal-to-noise ratios in the range S/N = 7- 
16. A radial velocity (RV) was obtained from each spectrum 
by using the c ross-correlation techn ique through a weighted 
binary mask (Baranne et al.||1996|. The mask is composed 
of more than 2000 sharp and isolated spectral lines in the 
CAFE wavelength range. The cross-correlation was per¬ 
formed in a ±30kms“^ range around the expected RV of 


the star. The peak of the cross-correlation function (CCF) 
was measured by fitting a four-term Caussian profile. This 
velocity was then corrected for the barycentric Earth radial 
velocity at mid-exposure time. Since we took several con¬ 
secutive spectra on the same nights, we decided to combine 
the RV values of the corresponding pairs in the cases where 
their individual signal-to-noise was low (i.e., S/N < 10) 
and mutual discrepancies were larger than 50ms“^. This 
procedure can also diminish the effect of short term vari¬ 
ability, such as the granulation noise, on the radial velocity 
(although here the expected amplitudes are of the order 
of few tens of ms“^). Some of the spectra were neglected 
due to low quality, mainly caused by weather conditions 
or problems with the stability of the ThAr lamp. We ob¬ 
tained 20 RV values in total which are reported, together 
with their observing times, in Table [l] and are compatible 
with the presence of a ~ 1 Mjup planet in the system (see 
Fig.and Sect.|^. 


Table 1: RV and BVS measurements of Kepler-539 from 
CAHA/CAFE. 


Date of observation 
BJD(TDB)-2450000 

RV 

(kms“^) 

errRv 

(kms“^) 

FWHM 

(kms“^) 

BVS 

(kms~^) 

6116.5291378 

10.408 

0.026 

-0.041 

9.696 

6124.3823828 

10.383 

0.029 

-0.065 

9.517 

6138.4327769 

10.387 

0.018 

-0.117 

9.656 

6523.4275256 

10.402 

0.024 

-0.100 

9.572 

6598.2879800 

10.448 

0.033 

0.081 

9.518 

6804.5209537 

10.439 

0.029 

-0.146 

9.250 

6811.4391381 

10.475 

0.035 

0.159 

9.407 

6821.6225006 

10.432 

0.049 

-0.118 

9.198 

6834.4194255 

10.437 

0.025 

0.086 

9.342 

6859.5818459 

10.409 

0.038 

-0.391 

10.036 

7136.6259467 

10.285 

0.039 

-0.266 

9.541 

7152.5198831 

10.431 

0.042 

0.087 

9.422 

7160.5096750 

10.450 

0.021 

-0.267 

9.312 

7169.5836757 

10.500 

0.032 

0.038 

9.153 

7193.4502198 

10.425 

0.045 

-0.088 

9.556 

7254.4473639 

10.361 

0.035 

-0.021 

9.727 

7258.3855966 

10.400 

0.045 

-0.106 

9.559 

7263.4322174 

10.402 

0.028 

0.014 

9.519 

7264.4220946 

10.413 

0.025 

0.018 

9.599 

7266.3737337 

10.449 

0.034 

0.169 

9.995 


2.3. Spectral analysis and the age of Kepler-539 

We derived the spectroscopic parameters of the host star 
Kepler-539 from the co-added CAFE spectra, which has a 
S/N ratio of about 40 p er pixel at 5500 A. Fo l lowing the 
procedures described in |Candolfi et al.| (|2013[ |2015|), we 
used a customised IDlj^ software suite to fit the composite 
CAFE spectrum to a grid of synthetic theoretical spectra. 
The latter were calculated with the stellar spectral s ynthe- 
sis program SPECTRUM ( Cray &: Corbally 199^ using 
ATLAS9 plane-parallel model atmospheres (Kurucz 1979), 
under the assumptions of local thermodynamic equilibrium 
(LTE) and so lar atomic abundances as given in Crevesse & 
Sauval (1998). We fitted spectral features that are sensitive 
to different photospheric parameters. Briefly, we used the 
wings of the Balmer lines to estimate the effective temper¬ 
ature Teff of the star, and the Mgi 5167, 5173, and 5184 A, 


^ The acronym IDL stands for Interactive Data Language and 
is a trademark of ITT Visual Information Solutions. 
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Fig. 1: The twelve transit events of Kepler-539 b observed by Kepler in long-cadence mode. The transits at epoch -5 and 
2 are clearly affected by star-spot-crossing events. Times are in BKJD (Barycentric Kepler Julian Date - equivalent to 
BJD(TDB) minus 2454833.0). 




Fig. 2: The three transit events of Kepler-539 b observed by Kepler in short-cadence mode. Times are in BKJD 
(Barycentric Kepler Julian Date - equivalent to BJD(TDB) minus 2454833.0). 


the Cai 6162 and 6439 A, and the Nai D lines to deter¬ 
mine the surface gravity logg*. The iron abundance [Fe/H] 
and microturbulent velo city Umicro was derived by a pplying 
the method described in Blackwell & Shallis (1979) on iso¬ 
lated Fei and Fell lines, lb determine the rnacroturbulent 
velocity Umacro, we adopted the calibr ation equations for 
solar like stars from Doyle et al. (2014). The projected ro¬ 


tational velocity v sin i* was measured by fitting the profile 
of several clean and unblended metal line^ We found that 
Kepler-539 has an effective temperature of Teff = 5820 ± 
80K, logg* = 4.4 ±0.1 (cgs), [Fe/H]= -0.01 ± 0.07dex, 
'ymicro = 1-1 i 0.1 kms“^, t’macro = 3.2 ± 0.6 kms“^, and 


^ Here i* refers to the inclination of the stellar rotation axis 
with respect to the line of sight. 
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Fig. 3: The entire Kepler light-curve data of Kepler-539. The large stellar variability can be reasonable interpreted as 
induced by a star-spot activity. Times are in BKJD (Barycentric Kepler Julian Date - equivalent to BJD(TDB) minus 
2454833.0). The red arrows mark the mid-times of the twelve transits of Kepler-539 b. 


V sini* = 4.4 ± 0.5 kms ^ (Table [^. According to the 


Straizys & Kuriliene (1981) calibration scale for dwarf stars, 
the effective temperature of Kepler-539 translates to a G2 V 
spectral type. 



X (A) 


Fig. 4: CAFE co-added spectrum of Kepler-539 (black line) 
encompassing the Lii 6707.8 A absorption doublet. The 
best fitting ATLAS9 spectrum is overplotted with a thick 
red line. The vertical dashed line marks the position of the 
Li doublet. The lowest part of the plot displays the residuals 
to the fit. 


The CAFE co-added spectrum of Kepler-539 reveals 
the presence of a moderate Lii 6707.8 A absorption dou¬ 
blet (Fig. |3. We estimated the photospheric lithium abun¬ 
dance of the star by fitting the Li doublet using ATLAS9 
LTE model atmospheres. We fixed the stellar parameters 
to the values given in Table and allowed our code to 
fit the lithium content. Adopting a correction for non-LTE 


effects of -1-0.006 dex ( |Lind et al.||2009| ), we measured a 
lithium abundance of A(Li) = log (n(Li)/n(H)) -|- 12 = 
2.48 ±0.12 dex. 

The photospheric lithium content and the moderate ro¬ 
tation period of Kepler-539 suggest that the star is rela¬ 
tively young. Fig.j^shows the lithium abundance of Kepler- 
539 compared to the A(Li) of the dwarf stars of the 
Pleiades, Hyades , and NCC752 open clust ers, a s listed in 
Soderblom et al. (1993), Pace et al. (2012), and Sestito et 
n" ( |2004p , respectively. Kepler-539 has to be older than 
~0.i Cyr, as the star falls below the envelope of the 
Pleiades. Although lithiu m depletion becomes in effective 
beyond an age of 1-2 Cyr (Sestito & Randich 2005), Kepler- 
539 lies between the envelopes of the other two clusters, sug¬ 
gesting an age intermediate between the age of the Hyades 
(-0.6 Cyr) and the age of NCC752 (-2.0 Cyr). This is 
further confirmed by the fact that the lithium content of 
Kepler-539 is intermediate between the average lithium 
abundance measured in early C-type stars of 0.6-Cyr-old 
open clusters (A(Li) = 2.58 ± 0.15 dex) and that of 2- 
Cyr-old open clusters (A(Li) = 2.33 ± 0.17dex; Sestito & 
Randich]|2005 ). 

We used Eq. (32) from [Barnes (2010) and the rotation 
period of Kepler-539 to infer its gyrochronological age, as- 
suming a convective turnover time-scale of Tc = 34 days 
(Barnes & Kim||2010) and a zero-age main sequence rota- 
tion period of Pq = 1-1 days (|Barnes||2010|. We found a 
gyrochronological age of 1.0 ± 0.3 Cyr, whicn supports the 
relatively young scenario. Our estimation is in good agree- 
ment with the 1.15 Cyr gyro chronological age predicted by 
Walkowicz & Basri ( 2013a|b|_and with th^ estimated using 
theoretical models (see Sect. I^and Table[^. 

By having short-cadence data over a quite large times- 
pan, it is in principle possible to perform an asteroseismic 
analysis and try to precisely estimate the age of a star. 
However, in the case of Kepler-539, the SC data are only 
available for less than half the quarters. In particular, the 
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continuous time baseline of SC data is not more than 125 
d in the best case (BJD2456015 — 2456139). Furthermore 
the activity variation is predominant in the data. We anal¬ 
ysed such SC data by filtering the variability due to activ¬ 
ity with moving averages, but we did not find evidence of 
a clear frequency comb in the power spectrum and there¬ 
fore we did not get particular clues about patterns due to 
solar-like oscillations. We concluded that no reliable aster- 
oseismic analysis could be performed for Kepler-539 based 
on the Kepler data. 


Table 2: Nearby visual companions around Kepler-539 
(from Adams et al. 2012) and their dilution effect on the 
depth of the transit events. 


Companion 

Distance (") 

Ks (mag) 

AK, 

Dilution 

B 

2.49 

23.2 

8.6 

0.005 % 

C 

3.56 

22.4 

8.0 

0.01 % 

D 

4.99 

22.7 

8.2 

0.008 % 

E"- 

5.94 

17.1 

4.0 

1.31 % 


Notes. ““Also known as KIC 6471028. 
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bution, see Figueira et al.|[20T6 1. This value indicates only 
a weak correlation between the two parameters, which is 
only significant at a 2cr level for the median value and still 
at less than 3tT for the upper boundary. Consequently, we 
assume no correlation between both parameters, although 
we warn about the weak probability of correlation. 


1 |— • Pleiades (Soderblom et al., 1993) 

■ Hyodes (Pace et al., 2012) 

A NGC 752 (Sestito et ol., 2004) 


3. Physical properties of the system 


6200 6000 5800 5600 

Effective temperature (K) 


5400 


Fig. 5: Lithium abundance s of F- and G-type dw arf stars 
in the Pleiades (black dots; Soderblom et al.|1993 ), Hyades 
(gre en squares; Pace et al.|2012 ), and NGC 752 (blue trian¬ 
gles; Sestito et al.||200j ) open clusters. The red dot marks 
the position of Kepler-539. 


2.4. Excluding false-positive scenarios 

A faint pulsating variable star or a binary system in the 
background/foreground can mimic a planetary transit sig¬ 
nal on the target star. High-resolution ima ges are extremely 
useful to exclude this possibility (e.g., Lillo-Box et al. 
2014c). A iCg-band, high-resolution, adaptive optics im- 


age of Kepler-539 was obtained w i th AR IES on the 6.5 m 
MMT telescope by Adams et al. (2012^ who foun d four 
stars within 6" of the target, see TaBIej^ Following Lillo- 


Box et al. (2012), we have estimated the dilution effect 
caused by each companion, finding that it is very small 
for the three closest targets, roughly 0.023% in total and 
thus negligible. Concerning the very bright companion “E”, 
this is another object targeted as KIC 6471028. As it lies 
outside the Kepler aperture, this star does not contami¬ 
nate Kepler-539. KIC 6471028 was also detected with the 
AstraLux North instrument mounted on the CAHA 2.2 m 


1.00 


0.99 


0.98 


0.97 



-0.003 -0.002 


-0.001 0.000 0.001 
Orbitol phose 


0.002 0.003 


telescope by Lillo-Box et al. (20121, who estimated that this 
is a K2-K4 background dwarf! 

An intense stellar activity could mimic the presence of 
a planetary body in the RV signal, thus causing a false pos¬ 
itive case. We have also analysed such a possibility by de¬ 
termining the bisector velocity span (BVS) from the CAFE 
spectra (the values are reported in Table[B. The bisector 
analysis provides a Pearson correlation coCTficient between 
the RV and the BVS of 0.28to;42 (median value and 95% 
confidence intervals from the posterior probability distri- 


Eig. 6: Phased Kepler long-cadence (top light curve) and 
short-cadence (bottom light curve) data zoomed around 
transit phase. The TTVs (see Sect.|^ were removed from 
the data before plotting. The jktebop best fits are shown 
using solid lines. The residuals of the fits are plotted at the 
base of the figure. 


For the determination of the physic al parameter s of th e 
Kepler-539 system, we proceeded as in |Ciceri et al. (2015). 
We first selected for analysis all data within two transit 
durations of a transit, ignoring the two with only partial 
coverage in the Kepler light curve, and converted from flux 
to magnitude units. Each transit was then detrended by 
fitting a polynomial versus time. High polynomial orders of 
3 to 5 were required to account for the brightness variations 
of the host star caused by spot activity. 

We then fitted the photometry and the CAFE RV 
measurements sim ultaneously using the jktebop code 
(Southworth 2013), after modifying it to allow fitting for 
individual times of mid-transit. The Kepler LC and SC 
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Table 3: Final parameters of the planetary system KOI-0372. 


Parameter 


Nomen. 


Unit 


Value“ 


Stellar parameters 

R.A. (J2000) . 

Dec. (J2000) . 

Spectral type** .... 


Kepler magnitude . 

.. 

.mag 

Effective temperature . 

. Teflf . 

.K . 

Iron abundance . 

....[Fe/H] .... 

.dex , 

Lithium abundance . 

.^(Li) . 

.dex , 

Microturb. velocity . 


... km s“ 

Macroturb. velocity'^ . 

.... 'Umacro .... 

... km s“ 

Proj. rotat. velocity . 

... .V sin u .... 

... km s“ 

Stellar rotation period'^ . 

. -Prot . 

.day , 

Stellar age (from gyrochronology) . 


.Gyr 

Stellar age (from isochrones) . 


.Gyr 

Stellar mass . 

.M* . 

.MfTl 

Stellar radius . 

. R* . 

. Rq 

Stellar mean density . 

. P* . 

. PQ ■ 

Stellar surface gravity . 

.log fir* . 

.cgs . 


Planetary parameters (Kepler-539 b) 

Planetary mass . 

Planetary radius . 

Planetary mean density. 

Planetary surface gravity . 

Planetary equilibrium temperature 
Safronov number . 


Orbital parameters 
Time of mid-transit .... 

Orbital period . 

Semi-major axis . 

Orbital inclination. 

Fractional star radius .. 
Fractional planet radius 

RV semi-amplitude . 

Barycentric RV . 

Eccentricity . 


19‘'56“29:39" 

41'’52'00.3" 

G2V 

12.39 

5820 ± 80 
-0.01 ±0.07 
2.5 ±0.1 
1.1±0.1 
3.2 ±0.6 
4.4 ±0.5 
11.769 ±0.016 
1.0 ±0.3 

1 1 + 0 . 6 + 0.3 
^•^- 0 . 0 - 0.1 
1.048 ±0.034 ±0.025 
0.952 ±0.017 ±0.007 
1.216 ±0.054 
4.502 ±0.014 ±0.003 


.Mb .... 

. Mjup .... 

0.97 ± 0.29 ± 0.02 

. Rh ■ ■ ■ ■ 

. . .... Pjup .... 

0.747 ± 0.016 ± 0.006 

. Pb .... 

. PJup . 

2.18 ± 0.66 ± 0.02 

. 9h ■■■■ 

. ms“^ _ 

43 ± 12 

T 

.K . 

387.6 ± 6.0 

. e . 


1.24 ± 0.37 ± 0.01 




■ To . 

Porb . . . . . 

.. a . 

. A . 

. Ta . 

..BJDtdb ■ • 

.days .... 

.an . 

.degree .... 

. Tb . 


■ Ka . 

.ms“^ .... 

■•7 . 

... kms“^ ... 

.. e . 



2 455 588.8710 ±0.0030 
125.63243 ± 0.00071 
0.4988 ± 0.0054 ± 0.0039 
89.845 ± 0.086 
0.01057 ±0.00078 
0.000854 ± 0.000085 
132.3 ±6.3 
9.959 ±0.007 
< 0 . 39 *^ 


Notes. Where there are two error bars, the first is a statistical error, coming from the measured spectroscopic and photometric 
parameters, while the second is a systematic error and is given only for those parameters which have a dependence on theoretical 
stellar models. 

“ The adopted parameters assume a circular orbit. 

^ With an accuracy of ± 1 sub-class. 


Using the calibration equatio ns of Doyle et al. (20141. 

From McQuillan et al. (20131. 

The 99% confidence upper limit on the eccentricity from a model in which the eccentricity is allowed to vary in the fit. 


data were fitted separately, accounting for the long effec¬ 
tive exposure times in the LC data by oversampling the 
fitted model by a factor of five. 

The following parameters were htted: the fractional 
radii of the two objects (r* = i?*/a and = Rp/a, where 
a is the orbital semi-major axis), the orbital inclination i, 
the time of midpoint of each transit, the velocity ampli¬ 
tude KT*, the systemic velocity Uy, and the coefficients of 
the polynomial for each transit. The orbital period, Porbj 
and reference transit midpoint, Tg, were fixed to the values 
measured from the transit times; these were used only for 
phasing the RV measurements. As in t he other papers of 
our series (e.g. Mancini et al. 2013a|b ), a quadratic limb 
darkening law was used, with the linear term fitted and the 


quadratic term fixed to 0.27 (Sing 2010). We rescaled the 
error bars of the LC, SC and RV data to give a reduced 
of Xy = 1-0 for each versus the fitt ed model. Based on ou r 
experience with similar studies (e.g. Mancini et al. 2014a|b |, 
this procedure is necessary for taking into account underes- 
timated data errors and additional sources of white noise; 
for example an RV jitter due to the activity of the host star. 


The best fits to the photometry and the RVs are shown 
in Figs. and 7^ The uncertainties in the fitted parameters 
were derived by running both Monte Carlo and residual- 
permutation simulations (Southworth 20081 and choosing 
the larger of the two error bars for each parameter. The 
results obtained using the SC data are more precise than 
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Fig. 7: Upper panel: phased RVs for Kepler-539 and the best 
fit from JKTEBOP. Lower panel: residuals of RVs versus best 
fit. 


those from the LC data, so were adopted as the final set of 
photometric parameters. 

We find that the observations are fully consistent with 
Kepler-539 b moving on a circular orbit. However, we per¬ 
formed the fit both fixing the orbital eccentricity e to zero 
and fitting for it via the combination terms e cos w and 
esinw, where w is the argument of periastron. In the latter 
case, we obtained e = O.lSto^ig- We then used the Bayesian 
information criterion (BIG) to evaluate the preferred sce¬ 
nario, finding that the eccentric case is not favoured over 
the zero eccentricity model. 

Finally, we estimated the age and the physical prop¬ 
erties of the system from the SC results and the T^g and 
[Fe/H] measured from our spectra. Constraints from the¬ 
oretical stellar models were used to make the solution de¬ 


terminate. We followed the HSTEP approach (Southworth 


2012 and references therein), which yielded measurements 
of the properties of the star and planet accompanied by 
both statistical and systematic errorbars. The final results 
are reported in Table|^ 


4. Planet Kepler-539 c from transit time variation 


We have analysed the 12 mid-transit times of Kepler-539 b 
to characterise the ephemeris of the transit and check if 
there is a possible transit time variation (TTV), which 
could be a sign of the presence of additional bodies in this 
system. Since the transits at cycles -5 and 2 are quite af- 
fected by starspots, we used the PRiSiv^and GEMc]^ codes 
(Tregloan-Reed et al. 2013 20151 that allowed us to fit both 
the fuff transit events and the snorter starspot-occultation 
events contemporaneously. The parameters of the starspots 
coming from the fits are reported in Appendix[A| The mid¬ 
transit times for the other epochs were estimated using jk- 
TEBOP. The resulting timings are tabulated in Table|^and 
were fitted with a straight line to obtain the orbital period, 
^orb = 125.63243 ± 0.00071 days, and the reference time 


® Planetary Retrospective Integrated Star-spot Model. 
^ Genetic Evolution Markov Chain. 


Table 4: Kepler times of transit midpoint of Kepler-539 b 
and their residuals. 


Time of minimum 
BJD(TDB)-2400000 

Cycle 

no. 

Residual 

(JD) 

54960.74947 ±0.00109 

-5 

0.040612 

55086.35000 ± 0.00040 

-4 

0.008712 

55212.96694 ±0.00042 

-3 

-0.006781 

55337.57977 ±0.00043 

-2 

-0.026380 

55463.23346 ±0.00118 

-1 

-0.005126 

55588.87808 ± 0.00041 

0 

0.007068 

55714.51794 ±0.00043 

1 

0.014495 

55840.13340 ±0.00096 

2 

-0.002478 

55965.77212 ± 0.00043 

3 

0.003803 

56091.39778 ±0.00030 

4 

-0.002963 

56217.00030 ±0.00030 

5 

-0.001233 

56342.66701 ± 0.00030 

6 

0.001403 


Table 5: The next mid-transit times of Kepler-539 b. 


Cycle 

BJD(TDB) 

Date 


Time (UT) 

16 

2457598.9899 

2016 Jul 

29 

11:40 

17 

2457724.6223 

2016 Dec 

02 

02:56 

18 

2457850.2547 

2017 Apr 

06 

18:07 

19 

2457975.8872 

2017 Aug 

10 

09:12 

20 

2458101.5196 

2017 Dec 

14 

00:29 


of mid-transit, Tq = 2 455 588.8710 ± 0.0030 BJD (TDB). A 
plot of the residuals around the fit (see Fig.l^ shows a clear 
and modulated deviation from the predicts transit times, 
meaning that the orbital period of Kepler-539 b is variable 
(the maximum TTV is 58.5 ± 1.6 min and was measured for 
the transit at cycle —5). The uncertainties of Porb and Tq 
have been increased to account for this, but we stress that 
these values should be used with caution as they are based 
on only 12 timings, of which two are incomplete and oth¬ 
ers are affected by starspot anomalies that can introduce 
offsets in the transit timing measurements - several studies 


have quantified this effect 

(e.g. 

Barros et al.||2013 |Oshagh| 

et al. 

2013 

Mazeh et al. 

2015 

loannidis et al. 2016|). We 


have reported the predicted mid-transit times m 'Tablej^for 
the next year and a half. 

In order to identify the origin of this modulated 
TTV, we used the TR ADES (TRAnsits a nd Dynamics of 
Exoplanetary Systems; Borsato et al.|20i4 | code, which al¬ 
lows modelling of the dynamics of multiple-planet systems 
by reproducing the observed mid-transit times and RVs, 
with the possibility to choose among four different algo¬ 
rithms. We thus investigated the possibility that Kepler- 
539 is a planetary system formed by a star and two plan¬ 
ets, b and c. We ran two series o f simulations, the fir st by 
selecting the PIKAIA algorithm ( Charbonneau| 1995), the 
second with the Pa rticle Swarm Optimization (PSO) algo¬ 
rithm ( |Tada|[2007| . In all the cases, we fixed tight bound¬ 
aries for the parameters of planet b (centred on the values 
in Table|^ and assumed the zero eccentricity case, which 
means that we fixed the argument of pericentre (wb) to 90°, 
and let only the mean anomaly (t'b) of planet b to freely 
float between 0 and 360 degrees. We use a time of refer¬ 
ence (or epoch, tepoch = 2455463.244) that is very close 
to the transit cycle 0, forcing the planet b to be initially 
at the pericentre. Because of this, all the solutions of both 
the algorithms returned zero values for uq. Eurthermore, 
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Cycle number 

Fig. 8: 0-C diagram for the timings of Kepler-539 b at mid-transit versus a linear ephemeris. The timings in blue refer 
to those coming from the Kepler LC data, while those marked with green are from SC data. The red points refer to the 
two LC transits affected by star spots. The two incomplete LC transits are marked with orange points. 


we assumed for simplicity the two planets to be coplanar. 
We investigated different configurations of two planet sys¬ 
tem, with the period of planet c. Pc, varying from 150 to 
2500 days. We fitted the pair (cc cos Wc; Cc sinwc) instead 
of (ccWc) to avoid correlation between the parameters. We 
set 0.9 as the highest value allowed for the eccentricity, 
while we let ujc and Vc freely vary as for i/b. The solutions 
obtained from the TRADES fitting processes were further 


analysed using the Frequency Map Analysis (FMA; Laskar 


et al 


19921 to check if they are stable, 
found that the best solutions are those related to 
configurations with Pc > 1000 days and a mass of the 
planet c between 1.2 and 3.6Mjup on a very eccentric or¬ 
bit (0.4 < Cc < 0.6). The final parameters, with confidence 
intervals estimated with a bootstrapping method, are re¬ 
ported in Table[^ for the five best solutions, based on the 
, and displayed in figures shown in Appendix[B| The es¬ 
timated confidence intervals are very tight to the fitted 
values because the algorithms found solutions that corre¬ 
spond to minimum surrounded by very high peaks in the 
X^-space. We also exhaustively investigated narrow region 
(200 — 800days) of Pc (in steps of about 300 days), but 
the solutions that we obtained are not favoured because 
the corresponding is much higher than those reported 
in Table[6l 

Knowing the orbital period of Kepler-539 c, we can eas¬ 
ily estimate the semi-major axis and the equilibrium tem¬ 
perature of the planet c for each solutions. These values 
are shown in right-hand columns of Ta ble[6l Moreover, we 
can use Eq. (7.29) from Haswell (|2010|) for estimating the 
maximum TTV signal expected for Kepler-539 b, i.e. 


TTVmax = (Ma + Mb)ec(ab/ac)^Pc- (1) 

The last column of Table[^ reports the values of TTVmax 
for each of the five solutions. The best TRADES solution is 
also that for which the maximum TTV signal is higher than 
the others and in better agreement with the observations. 

Very interestingly, the high eccentricity of planet c gives 
a direct explanation of the modulation of the TTV of planet 
b: when planet c moves close to planet b, it gravitationally 
kicks its smaller sibling, causing the maximum in the TTV 
signal; after the conjunction, planet b re-circularizes own 
orbit very slowly, because of tidal interactions with the star 
and planet c moving far from it. This causes the decreasing 
of amplitude of its TTV that we see in Fig|^ 

We again stress that all the values presented in Table[^ 
should be taken with caution, because the existence of 


planet c can be definitively constrained only by new tran¬ 
sit observations, able to completely cover the TTV phase 
of planet b. In particular, one should verify if there is a 
periodic repetition of the peak and the modulation. 


5. Discussion and conclusions 

Thanks to precise RV measurements obtained with the 
high-resolution spectrograph CAFE, we confirmed the 
planetary nature of Kepler-539 b, a dense Jupiter-like 
planet with a mass of 0.97 ± 0.29 Mjup and a radius of 
0.747 ± 0.018 Rjup, revolving with a period of 125.6 days 
on a circular orbit around a G2 V star, similar to the Sun. 
The parameters of the parent star and the planet b were ob¬ 
tained by analysing the CAFE spectra and through a joint 
fit to the RV data and Kepler transit light curves. Orbiting 
at ~ 0.5 au from its host, Kepler-539 b is located not so 
far (« 0.45au) from its habitable zone (HZ), which has a 
width of ~ 0.72 au, as we estimated it based on the stel- 


by Kopparapu et al. 

( 

20131 |2014D. The parent star is quite 

active, as is shown 

oy the 0.047 mag peak-to-peak modu- 


lation present in the long time-series photometry and from 
the star-spot anomalies clearly visible in two of the transit 
events monitored by Kepler. Even though it has physical 
characteristics that resemble those of the Sun, Kepler-539 
is much younger: 1.0±0.3 Gyr as estimated from its lithium 
abundance and gyrochronology. 

An analysis of the Kepler mid-transit times of Kepler- 
539 b revealed a clear modulated TTV signal. We studied 
this with TRADES, a code able to make a simultaneous 
fit of RV and mid-transit times and compare the results 
with simulated data of various multi-planetary configura¬ 
tions. We found that the TTV can be explained through 
the presence of a 1.2-3.6Mjup Jovian-like planet c on a 
very eccentric (e = 0.43-0.61) and wider orbit than that of 
Kepler-539 b, with a period larger than 1000 days. Each of 
the five solutions that we have presented are related to a 
planetary system which is stable according to the results of 
an FMA analysis. 

Since the orbit of Kepler-539 c is highly eccentric, then 
the modulation of the TTV signal of Kepler-539 b can be 
easily explained. As the distance of planet c from the parent 
star varies with time, there is a periodic change of the mid¬ 
transit time of planet b. In particular, when the planet c is 
in conjunction with planet b, the gravitational interaction 
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Table 6: Left-hand columns: parameters of Kepler-539 c from the best fit of the Kepler mid-transit times and CAFE-RV 
measurements. The best five solutions obtained with TRADES are shown. Left-hand columns: parameters determined 
from the previous. 


Algortihm Me (Mjup) Pc (day) 




<Vc(°) 


a (au) 


Teg (K) TTVmax (min) 


PIKAIA 

O.DU_o 23 

1040t?9 

33.95l°;^« 

0 

u.^oz_o 006 

270.671°;°® 

34.2 

2.04 ± 0.04 

271 ±4 

31.0 ± 3.6 

PSO 

-I orr + 0.05 
-^•^'-0.05 

1705tl° 

90 1 q+0.84 

zo. J.o_o 72 

U.DUO_o 005 

166.27l°-|° 

39.5 

2.84 ± 0.08 

230 ±4 

9.3± 1.1 

PIKAIA 

-L-^U_o 11 

9451" 

I8.00t°;” 

u.4^o_0_oio 

329.2712 54 

57.2 

1.91 ± 0.02 

280 ±3 

11.8± 1.4 

PIKAIA 

9 79+0-15 
'^-0.19 

1784t|6 

349.141°;®° 

U.OD9_o 010 

296.56flll 

62.9 

2.92 ±0.08 

226 ±4 

17.8 ± 2.2 

PSO 

C) QQ+0. 15 


355.50l°;°| 

U.‘4O<'0.009 

343.541);®^ 

66.6 

1.94 ± 0.05 

278 ±4 

22.4 ± 2.7 


Notes. We used two different algorithm: the PIKAIA algorithm (|Charbonneau[|1995[) and the Particle Swarm Optimization (PSO) 


algorithm (Tada 
processes, while 


20071, respectively. The first five columns refer to the values of the parameters determined from the fitting 
he others values have been determined using the previous ones. In particular, St represents the maximum TTV 


for the inner planet b. 



1 10 100 1000 

Period (days) 


Fig. 9: The period-planetary mass diagram for transiting planets in the mass range O.SMjup < Mp < llMjup. The size 
of each circle is proportional to the corresponding planetary radius, while the colour indicates equilibrium temperature 
(data taken from TEPCat). The dashed line demarcates the high-mass regime {Mp > 2Mjup). The positions of Kepler- 
539 b and Kepler-539 c are highlighted. Kepler-539 c is marked with a box since we do not know its radius. The error 
bars have been suppressed for clarity. 


between the two planets is maximum and we see a large 
amplitude of the TTV signal. Then, when planet c moves 
away, planet b re-circularizes its orbit very slowly and we 
see the gradual decreasing of the amplitude of its TTV. 

Given its relatively young age, the physical parame¬ 
ters of the Kepler-539 planetary system can be very pre¬ 
cious for astrophysicists working on theories of planet for¬ 
mation and evolution. If we visualise the known transit¬ 
ing planets, reported in TEPCal]^ in a period-mass dia¬ 
gram, we can see that Kepler-539 b and Kepler-539 c occupy 
sparsely-populated regions of this plot (here for the mass 
of Kepler-539 c we have considered the best value found by 
TRADES). 

We also underline that the transits of Kepler-539 b are 
quite deep (0.8%) and that the host star is relatively bright 
{V = 12.5 mag) and therefore amenable for exoplanet at¬ 
mosphere studies. The star-planet distance is large enough 
to consider Kepler-539 b as a non-inflated planet with a low 
irradiation level. Since there are only a handful of transiting 


® TEPcat (Transiting Extrasolar Planet Catalogue) is ava il- 
able at www.astro.keele.ac.uk/jkt/tepcat/ ( Southworth|20 iT I. 


Jupiter-like planets with a low irradiation level, this plane¬ 
tary system could be a precious target for future studies of 
exoplanet atmospheres of normal giant gas planets. 

Unfortunately, the parameters of Kepler-539 c are not 
well constrained as they are based on a poorly-sampled 
TTV (only 12 transit timings). More transits of Kepler- 
539 b and more precise RV observations with a larger- 
aperture telescope are required to properly model its TTV 
signal and accurately determine the existence and the phys¬ 
ical properties of Kepler-539 c on a quite solid ground. 
While there are now several high-resolution spectrographs 
with RV performance better than CAPE, the observations 
of new transits of Kepler-539 b are not straightforward with 
ground-based facilities. The 125.6-day orbital period and 
the long transit duration (9.31 hr) put serious limitations 
for the successful of a such long-term monitoring and re¬ 
quest a large observational effort with various telescopes 
located at different Earth’s longitudes. Furthermore, the 
use of different telescopes could compromise the precision 
of the measurements of Tq at the end. On the other hand, 
the observations of new mid-transit times of Kepler-539 b 
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can be easily performed by the incoming space telescope 
CHEOPS. This facility is therefore highly recommended for 
a new detailed follow-up study of the Kepler-539 planetary 
system. 
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Appendix A: Starspot parameters 

In this appendix we report a Table containing the starspot 
parameters that were determined from the prism+GEMC 
fitting of the transit light curves at cycles —5 and 2, as 
defined in Sect. 12.II 
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Table A.l: Starspot parameters derived from the prism-|-GEMC fitting of the transit light curves at cycles —5 and 2. 
“The longitude of the centre of the spot is defined to be 0° at the centre of the stellar disc and can vary from —90° to 
90°. ^The co-latitude of the centre of the spot is defined to be 0° at the north pole and 180° at the south pole. “Angular 
radius of the starspot (note that an angular radius of 90° covers half of stellar surface). ‘^Spot contrast (note that 1.0 
equals the brightness of the surrounding photosphere). 


Cycle 

6l(°)“ 


r.pot(°)° 

Pspot 

-5 

-26.18 ± 3.73 

28.31 ± 10.14 

36.43 ±6.80 

0.69 ±0.11 

2 

18.88 ± 11.79 

34.27 ± 8.21 

18.60 ±6.46 

0.70 ±0.28 
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Appendix B: TRADES simulations 


In this appendix we report plots based on the results of the 
TRADES simulations for modelling the TTV detected in 
the transit timings of KOT327b. In particular, the figures 
show RV plots (bottom panels) and 0-C diagrams from lin¬ 
ear ephemeris for planet Kepler-539 b (top panels). In both 
cases, the observed data are compared with those obtained 
from the simulations (see Sect.[^ The figures are ordered 
as in Table|^ Figures 1 , 3 and 4 were obta ined by using 


the PIKAIA algorithm (Charbonneau 1995), while figures 
2 and 5 with the P article Swarm Optiinization (PSO) algo¬ 
rithm ( Tada|[2007 ). 


I observations ° simulations 
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Fig. B.l: 0-C diagram (with residuals) from linear 
ephemeris for planet Kepler-539 b (top panel); observations 
plotted as black points (with error bars), simulations plot¬ 
ted as open blue circles. Bottom panel shows the RV ob¬ 
servations as solid black circles, simulations at the same 
BJD(UTC) as open blue circles, and the dotted gray line 
is the RV model for the whole simulation. This the best 
solution found by trades. 
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Fig. B.2: Same as for Fig |B.l| but showing the 2“^^ best 
lution found by trades. 
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Fig. B.3: Same as for FiglBT] but showing the 3’''^ best so¬ 
lution found by trades. 
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Fig. B.4: Same as for Fig |B.l[ but showing the 4*^ best so¬ 
lution found by trades. 
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Fig. B.5: Same as for Fig 
lution found by trades. 
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